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Abstract 

In the landscape, if there is to be any prospect of scientific prediction, it is crucial that 
there be states which are distinguished in some way. The obvious candidates are states 
which exhibit symmetries. Here we focus on states which exhibit discrete symmetries. Such 
states are rare, but one can speculate that they are cosmological attractors. We investigate 
the problem in model landscapes and cosmologies which capture some of the features of 
candidate flux landscapes. In non-supersymmetric theories we find no evidence that such 
states might be cosmologically favored. In supersymmetric theories, simple arguments sug- 
gest that states which exhibit R symmetries might be. Our considerations lead us to raise 
questions about some popular models of eternal inflation. 



1 Introduction: Distinguished States in the Landscape 



For almost 30 years, inspired by string theory, particle theorists have imagined that there is a 
sensible program of attempting to infer the laws of nature from some grand set of principles — 
from the top down. The emergence of the string landscape, however, has thrown this paradigm 
into question. If there is any truth to the idea of a cosmic landscapc[l, 2, 3, 4, 5], there is no 
realistic program in which one finds the string ground state which describes the world around 
us. At best, one can hope to do statistics, looking for strong correlations between known facts 
and laws of nature which can be studied in future experiments. 

At present, even an enthusiast must acknowledge that the landscape is more of a vision than 
a well-defined theoretical structure. At best, we have explored only a small subset of the possible 
set of states [6, 7]; most of these are supersymmetric or approximately so. The most interesting 
states, with broken supersymmetry and small cosmological constant, can hardly be said to have 
been understood in a systematic way. In the original KKLT construction, for example, while 
there are small numbers, there is not a formal, small parameter in which one can develop a 
systematic approximation scheme. On the other hand, the KKLT analysis seems plausible, 
and the statistics which arise from this and similar constructions seem robust. For example, 
the statistics discussed by Douglas and Denef[8, 9] emerge from simple assumptions about the 
low energy degrees of freedom, as well as assumptions of uniform distributions of lagrangian 
parameters with respect to some natural measure [10]. So it is perhaps not unreasonable to 
think that these model landscapes may provide reliable insights into the larger structure. 

It is not at all clear whether the (nearly) supersymmetric states are representative; serious 
arguments have been advanced that they are not[9, 11]. So, even accepting the validity of the 
landscape framework, it is unlikely that we will soon be able to reliably enumerate all possible 
classes of landscape vacua. Thus if there is to be any hope of prediction, it must emerge from 
some more general considerations. Among classes of states which have been studied, those 
which might be distinguished are: 

1. (approximately) supersymmetric states 

2. states with large internal volume 

3. states exhibiting warping 

4. states exhibiting (discrete) symmetries 
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While most thinking about states of string theory assumes some approximate supersymnie- 
try, it seems reasonable to think that supersymmetry is special, and that there might be vastly 
more supersymmctric states than non-supersymmetric statcs[ll, 9, 12, 13]. One can ask, on the 
other hand, whether cosmologically, such states might be favored. In [15], it was pointed out 
that primitive cosmological considerations might favor approximately supersymmctric states. 
Typical stationary points of an effective action are extremely unstable - they are not states 
at all. Approximately supersymmctric states, however, are inevitably highly metastable. Of 
the various distinguishing characteristics listed above, apart from supersymmetry, only large 
volume generically provides stability, and, at least within JIB landscapes, large volume does 
not appear generic. We won't say more about large volume or warping in this paper, as our 
considerations don't seem to provide any arguments that these states are special. 

In speculating about physics beyond the Standard Model, it has been a widely held belief 
that discrete symmetries might plausibly - and naturally - solve a range of problems. Such 
symmetries often arise at special points in the moduli spaces of critical string constructions. 
But in the framework of the flux landscape, discrete symmetries would seem to be rare[14]. 
The problem is that the vast number of states purported to exist within the landscape arise 
because of the existence of a large number of fluxes which can take a large number of values. 
However, in order to obtain symmetries, many fluxes must be set to zero, greatly reducing the 
"dimensionality" of the space. In other words, only an exponentially small fraction of states 
exhibit symmetries. 

Again, we can ask whether cosmology might in some way favor (or not) states exhibiting 
discrete symmetries. Here, the issue is not one of stability, but rather: are the symmetric 
states, within some realization of eternal inflation, often attractors. We might mean by this 
that the measure in some inflationary model is highly peaked at the symmetric states, but we 
will shortly formulate the question in a more naive and primitive way. 

Given that there is so much we do not understand about the landscape, the burden of 
this paper is not to prove that symmetric states are attractors, but rather to establish that 
this possibility is, or is not, plausible. For this, we seek to adopt models of the states within 
the landscape, the transitions between them, and the mechanisms by which they are popu- 
lated. For the states, the challenge is to study the highly non-supersymmetric states which 
are the antecedents of the would-be states of interest. We will consider a number of models 
for the states, starting with BP, and then continuing to "states" found as stationary points 
of supergravity models. These models will incorporate some, but not all, of the features we 
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might expect to hold in "reahstic" theories. In eonsidering transitions, as we have indicated, 
we win suppose radii and couplings to be such that these states are states, but that they are 
highly metastablc, so that they do not experience periods of exponential growth. One of our 
more drastic assumptions will be to treat the fluxes as quasi-continuous. In this way, we can 
think of tunneling as giving rise to a flow within the landscape, and ask whether, for a given 
starting point, the system flows towards symmetric states. Such a continuous flux approxima- 
tion has proven useful in thinking about landscape statistics [16], and has been applied in other 
cosmological contexts[17]. While surely unrealistic for this problem, it is hard to believe that 
if symmetries are not attractors in this framework that they are in any realistic setting. 

Finally, we need a model for cosmology. Ultimately, we must face up to questions of 
measures for eternal inflation. But our goal, here, is simply to establish whether there is any 
reason to think that the rare states in a landscape exhibiting symmetries are somehow favored. 
For this, again, the burden of this paper is not to establish conclusively that this is the case 
in an underlying, complete theory of gravity, but simply to establish some general conditions 
under which symmetries might plausibly be favored. Postulating a landscape with a large 
number of (very) metastable de Sitter states, we will consider possible "initial conditions" in 
which the universe starts in one such state. Then we can ask whether, for a non-negligible 
subset of possible starting points, the system finds its way to the symmetric state. To model 
this problem, we will suppose that the antecedents of the symmetric state are short lived and 
do not experience long periods of exponential growth. 

In the next section, we lay out certain basic issues involving symmetries in such a picture. 
We introduce the notion of a neighborhood in the landscape, and argue that it is plausible that 
an order one fraction of states lie in the neighborhood of states exhibiting discrete symmetries. 
In section three, we present simple arguments for our basic results: ordinary symmetries are 
not attractors, but R symmetries might be. We survey a set of model landscapes in section 
four. In this discussion, we enumerate some of the present obstacles to developing a landscape 
cosmology. We discuss the problem of symmetric points within model landscapes. The T^jZ^ 
orientifold is introduced as a model where the symmetry structure is readily understood. Our 
studies reinforce the arguments of section two, that points exhibiting R symmetries can readily 
be attractors, but other types of symmetries are not singled out. 

In section five, we discuss our conclusions. We explain that these studies are, at best, 
preliminary and suggestive. Making definitive statements about cosmology certainly requires 
a deeper understanding of non-supersymmetric configurations within the landscape. But our 
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discussion, we believe, makes clear that the viewpoint, advanced by some anthropic advocates 
(e.g. [18]), that symmetries play no role, docs not necessarily follow from adopting a landscape 
framework. There are clear ways in which supersymmetry and discrete symmetries might well 
emerge from cosmological considerations. Obviously continuous gauge symmetries are of great 
interest; we leave their study for future work. 



2 Neighborhoods of Symmetric Points 



States exhibiting discrete symmetries are rare in the flux landscape. But the fraction of states 
lying nearby such symmetric points need not be small. To understand this, suppose we have 
some number of symmetric fluxes, Ni, i = and symmetry-violating fluxes, Ua, a = 

1 . . .A. We want to ask whether the set of states in the neighborhood of the symmetric points, 
i.e., in which n < N, is an order one fraction of the total number of states - or more precisely, 
initial conditions. To model this, let's suppose that the initial state is taken to lie on the 
"hypersphere" 

Nf+nl = R' (1) 
(with summation oni,a) and ask what fraction of the volume of the sphere satisfies 

a'Ni > nl (2) 



For simplicity, we'll take the measure in the space of fluxes to be flat. So we need to 
compare: 

j d^Nd^nSiVN'^ + V? - R^) (3) 

with the same integral, restricted by ■n? < a^N"^ (the constant a, should be less than one, but 
our working assumption is that it need not be much less than one). To do the integrals, we need 

d 

to recall that the solid angle in d dimensions is: Cld = rp^- integral, in the unrestricted 

case, is simply 

r(^). 

With the restriction on the smallness of n, the number of states is 

J d^N J d^nSiVN^ + R^)e{a^\N'^\ - \n'^\). (5) 
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r(A/2)r(//2) 



Jo 



In the case that a ^ 1 we can readily estimate the integral, approximating the (1 — x)^/^ ^ 
factor in the integrand by 1. Then the integral behaves as a"^, while the total number of states 
behaves roughly as ^ ! so the overall suppression is a"^ ^4) ^ • ^o'^ large A and small a, this 

can be very substantial. If a ~ 1, on the other hand, numerical study of the integral indicates 
that the suppression with increasing A is quite modest. At the moment, the correct value of a 
is something we can only guess, but we will assume, in what follows, that it is not too small. 

One interesting question is how much disfavored are larger symmetries.^ As an example, 
wc consider the case, numerically, where we take a = 1, and compare I = 2,A = 10, with 
I = 6,A = 6. The difference in the two cases is a factor of 16, i.e. not huge. So it is not 
clear that there is an overwhelming preference for smaller symmetries (or at least for small vs. 
slightly larger symmetries). 



3 Symmetric Points: Attractors Within a Landscape? 

Having established that an order one fraction of states might lie in the neighborhood of a 
symmetric point, we ask: with our assumptions about initial conditions, is the system more 
likely to end up in a symmetric state than a non-symmetric state? 

Consider, first, the general question of tunneling within a flux landscape. In [15] this 
question was discussed from a slightly different point of view, which we will adopt here. The 
general idea is that the large number of states in the flux landscape arise because there arc a 
large number of possible fluxes, Ni = 1, . . .b, with each of the Ni taking on a large number of 
possible values. The potential, V{(p), is typically quadratic in fluxes (e.g. in the IIB landscape), 
V ~ N^, so: 



1. The difference in energy between states differing by a single unit of flux is of order A'^ 

2. The tension of the branes separating states differing by a single unit of flux does not 
depend significantly on the total flux. This follows from considering wrapped branes 
(NS5, D5). Alternatively, from the point of view of the fields in the low energy theory, 

^We thank L. Susskind for stressing the significance of this issue. 
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the tension is of order y^V{(f))A(l) ~ 1; here Acp, the change in typical fields (e.g. complex 
structure moduli in the IIB case) is of order 1/N. 

3. The bounce action for transitions between neighboring states is typically small for large 
N, scaling as 

S,{AN = 1)^^. (6) 

Ref. [15] focussed on the question: if a state happens to have small cosmological constant, 
is it stable against decay to its numerous neighbors with negative cosmological constant. Here 
we are asking a different sort of question. We are studying the ancestors of the symmetric 
states. Since we are interested in the symmetric states as candidate theories of nature, we will 
assume that they have small cosmological constant; the ancestors have positive cosmological 
constant. The estimate above indicates that these states are typically highly unstable - they are 
not states at all. But, as a model, we will take this type of formula literally (one can suppose, 
for example, that the states, for some reason, all have moderately large volume, enhancing the 
bounce action). In the same spirit, we will treat the fluxes as large enough that they can be 
treated as quasi-continuous [16]. In such a picture, there is a well-defined notion of neighborhood 
in the landscape. 

As a further, rather drastic simplification, we will ignore gravity in estimating decay am- 
plitudes. This is often valid (as discussed, for example, in [2]), for some range of fluxes and 
values of moduli, but it is not likely to hold generally. For example, among the many numerous 
metastable dS states, it is likely that some are long-lived, giving rise to (eternal) inflation. As 
we will see, even with such drastic assumptions, only in a limited set of instances do symmetric 
states seem singled out by cosmological considerations. 

3.1 Non- Symmetric States 

Suppose that the putative state of interest, with small cosmological constant, lies at some point 
in the flux lattice Nq. In this case, in the continuous flux approximation, we would expect the 
energy, near that point, to have the form: 

Eo = CaAiV« (7) 

Recall the basic assumption that there is a large, multi-dimensional space of states. So, from 
the tunneling perspective, with our assumption about initial conditions, it is clear that from 
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most starting points, one will "miss" the state Nq, ending up eventually tunneling to negative 
cosmological constant big crunches. The chance of ending up in the "correct vacuum" will be 
of order one divided by the number of low cosmological constant states. 

3.2 The neighborhood of symmetric states 

Consider a state which is symmetric or approximately symmetric under an ordinary discrete 
symmetry. Some subset of fluxes, Ni, i = 1,...,B, respect the symmetry (they are neutral 
under the symmetry, and there are minima of the resulting potential in which only fields neutral 
under the symmetry have expectation values). The rest, n^, a = 1, . . . , A, break the symmetry. 
The putative, low cosmological constant, symmetric state we denote, again, by Nq. The set of 
states with flux 

\na\<\Ni-Nf\ (8) 

define the neighborhood of the symmetric point. We will take n to denote an average Ua, N slu 
average N^. 

If we assume that in the symmetric state there are no moduli, then in the continuous flux 
approximation, the energy landscape nearby has the form 

/ A 
Eo= fiJ^i^J + E dabriarib (9) 

i,J=l a,b=l 

Again, in the continuous flux approximation, we can diagonalize g. We can immediately dis- 
tinguish two cases. 

1. gab has only positive eigenvalues. In this case, if one starts in the neighborhood of the 
symmetric state, at least the transitions which change the Ua's will tend towards the 
symmetric state. A priori, we might simply assign a probability 1/2 for the possibility 
that a particular eigenvalue is negative, so there is (compared to factors like 10^°°) a 
modest suppression of this possibility. 

2. gab has some negative eigenvalues. In this case, the corresponding n's will tend to grow, 
and the system will not tend towards the symmetric state. 

But even in the case where all of the eigenvalues of g are positive, there is nothing partic- 
ularly special about transitions to the symmetric Ha = state. Assuming the eigenvalues of g 
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are all similar, if na Ni, then 



Sb{5na = 1) 
Sb{6Ni = 1) ~ n3 



(10) 



So transitions which change n are much slower than those which change N. Under such cir- 
cumstances, the notion of a neighborhood is not particularly relevant to the tunneling process. 
Unless one starts in the state with Ua = 0, one will not reach the symmetric states; one may 



Suppose, on the other hand, that the elements of g, in addition to being positive, were 
far larger that those of /. Then transitions involving changes in n would be faster than those 
involving changes in N, and states in the neighborhood would, preferentially, end up in the 
symmetric state. Such a situation could conceivably arise if there were approximate moduli 
in the symmetric state. Then, in the nearby states, some quantities (volumes of cycles, for 
example) might be quite different than in the symmetric state, and tunneling rates could be 
quite different. We will a phenomenon of this sort in the case of R symmetric vacua. 

3.3 R Symmetries 

Let's turn our attention to R symmetric states. Here there is an important difference with the 
non-R case. Changes in N, in the symmetric limit, are not associated with changes in energy, 
so they are potentially highly suppressed. 

To go further, we need to think more carefully about the R symmetric states. In general, 
suppose that we have fields, Xj, which transform under the R symmetry like the superpotcntial, 
and fields, (j)a, which are neutral. There are additional fields, pi, which transform diflFerently 
than the X's under the R symmetry. Then the superpotcntial has the form 



There will be supersymmetric, R symmetric solutions (X = p = 0) provided there are more 4> 
type fields than X type fields. In this case, in general, there will be some number of moduli. 
This condition holds in many Calabi-Yau compactifications, for example[14]. 

In the present context, this means that classically the vacua with Ua = possess moduli. 
Our working assumption will be that these moduli are fixed by non-perturbative dynamics, as 
discussed in [19], with supersymmetry broken at comparatively low scales. For non-zero Ua, the 
situation is quite different. We can study an effective action for the moduli. Their potential is 



reach a state (iVo,na), but then one will transition to big crunches with negative cc. 



W = XiFi{cl)o,,Pi) 



(11) 
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of order n^. Correspondingly, the masses-squared of the hght moduh are of order ri^. As we 
will see in model examples shortly, the corresponding ground states can be dS or AdS. If much 
of the neighborhood of the symmetric point is AdS, we don't expect a significant probability 
to reach the symmetric point. So we will assume that the neighbors of the symmetric point are 
predominantly dS. 

In order to compare tunneling rates for processes with changes in N with those with 
changes in n, we need to understand how the properties of the lagrangian for the light moduli 
depend on N. We expect that this dependence is quite weak; indeed, we might guess that it 
goes as n^/N"^. We can illustrate the point, again, with a simple field theory model. Consider 
a theory with three fields, X and 0, %, where X transforms like W and 0, x, are neutral. For 
the superpotential we take 

W = XX{(i> + x){4>-x) (12) 

(the precise dependence on <p and x is not crucial, but makes the equations simple). There is 
a moduli space, on one branch of which 

(f) = X = v (13) 

On this branch, there are two massive fields, X and $ = -^((p — x) , and one massless field, 
-0 = '^{4> + x)- Now, the effect of some underlying theory in which we turn on fluxes, Ui, would 
be to generate a superpotential of the form, say 

AVF = e(a + /?^ + 7$). (14) 

where e ^ n. Such a potential can fix the remaining modulus (and break supcrsymmctry). It 
also generates a tadpole for of order e^, and a A-depcndcnt term in the energy of order 

We will see some of these features in more detailed models in subsequent sections. 

We will be interested in the relative size of transition rates for changes in n and changes 
in N. We have understood the energy splittings between the different types of states; we need 
to ask about the tensions. One might think that the tensions for changes in N would also 
be suppressed, as, for BPS configurations, the tension is proportional to /SW . But already in 
field theory, domain walls separating R symmetric states are not BPS. As a simple example, 
consider a theory with superpotential 

W = X{(j)'^ - ij?) 
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This model has an R symmetry and a Z2 symmetry. The states (p = ^fJ- are separated by a 
region in which (p ^ 0, and 7^ 0, but in which d^X = 0, so they are not BPS, and their 
tension is certainly non-vanishing. In the flux case, we expect, again, that the tension is not 
very sensitive to N. So the bounce action for changes in A'^ is larger than that for changes in n 

by 

^^(^^ = 1) = ^, (16) 
SbiAn = 1) n9 ■ ^ ' 

So in this case, transitions towards the symmetric point are likely to be much faster than other 
transitions. 

This discussion, while suggestive, is not complete. Much depends on the states in the 
low energy theory, defined by the n^'s. For example, this theory may itself have many AdS 
states, as we will see, which can mean that the system does not reach the symmetric point. 
We will discuss this issue in the models of the subsequent sections. Still, we see a sense in 
which symmetric points, and especially R symmetric points, are special, and could plausibly be 
attractors in a flux landscape. 



4 Model Landscapes 



We are far from understanding the structure of any realistic landscape. We will review some 
of the reasons for this shortly, but just note for now that this means that any understanding 
we have of the landscape - especially in a cosmological setting - comes from consideration of 
model landscapes, which it is hoped capture some of the important features of what might be 
the "real" situation. 



4.1 The Bousso-Polchinski Model 



The first compelling model of a landscape was that of Bousso and Polchinski (BP) [2]. This 
was based on consideration of four-form fluxes in Type II theories. One notable feature of the 
model is the assumption that all moduli are flxed for any choice of fluxes, in every state (and 
flxed, essentially to the same value). The energy spectrum is taken to have the form: 

Eo = ^Niqf-Ao. (17) 

The Qi's are constants, independent of TVj. They are assumed to be small enough that all 
tunneling amplitudes are small. This requires (in accord with the estimates in [15]) that the 
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internal manifold be large, with volume scaling as a positive power of the flux. This model is 
extremely useful, first, for illustrating the idea of a diseretuum: the model exhibits a nearly 
continuous distribution of energies for large fluxes. It also provides a model for eternal inflation. 

With these assumptions, if the q^s are comparable, transitions tend to decrease the fluxes 
uniformly. One can make a rough estimate of the transition rate by ignoring gravity and using 
the formulas appropriate to the thin wall limit[20, 21]. This was done in [2], in the limit that 
the compactification volume and the radii of various three cycles are large. The spirit of our 
analysis is to consider the case that radii are of order the fundamental scale, and to examine 
scaling with N. Then the bubble wall tensions are of order 1, while the energy splitting between 
states is of order N. So the bounce action behaves as As noted in [2], even moderately 

large volume can significantly suppress the tunneling rate. Note that if one of the NiS becomes 
much smaller than the others, further transitions in that flux are suppressed, i.e. all fluxes tend 
to decrease uniformly. 

4.2 IIB Compactifications: Limitations of the BP Model and The Continu- 
ous Flux Approximation 

Known flux models differ from that of BP in some significant ways. First, when moduli are 
fixed, they are typically non-trivial functions of the fluxes (this, and related points, have been 
made in [17]). In the IIB theory, for example, the potential for the moduli assumes the usual 
N = 1 supergravity form: 



where i,j are indices labeling the dilaton, complex structure and Kahler moduli. The superpo- 
tential, W , is linear in the fluxes, so 



where Ni,Nj are the fluxes and (j)a denote the various fields. If we can think of all fluxes as 
large, it is natural to write: 



Then, at stationary points of the potential, the (f)a are independent of N, and the potential is 
proportional to N^. But both (pa, and Vq, the values of the potential at the stationary points, 
are dependent on the Ui's. We will illustrate this with some simple toy models shortly. 




(18) 



(19) 



Ni = Nui. 



(20) 
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In the IIB models, there is another important difference: there are no flux-independent 
terms in the classical potential derived from the supergravity lagrangian. As a result, there is no 
obvious identification of the parameter or positivity of the parameters q^, and no clear sense 
in which the fluxes should tend uniformly towards zero with successive tunneling transitions. 

But perhaps the most crucial difference is related to the Kahler moduli (which we will de- 
note generally by p). The potential for the complex structure moduli has stationary points, but 
at large radius of the internal manifold, there are no stationary points for p; both perturbatively 
and non-perturbatively, the potential tends to zero for large p. As a result, the system is prone 
to run off to large radius. The value of p, and of the potential (and other quantities) is not 
similar at the KKLT point and at adjacent points. There may be stationary points (metastable 
vacua of the full quantum system) at small p (i.e. p ~ 1) for some or most choices of flux (or 
slightly larger p, as in KKLT, see below), but one would expect that if one starts in a high 
cc state with all moduli fixed, one is likely to tunnel to a configuration which simply rolls out 
to large p. Understanding the significance of this (assuming it has a sensible interpretation) 
probably requires a deeper understanding of eternal inflation, and especially how states are 
populated. 

4.3 KKLT and Limitations of the Continuous Flux Approximation 

The KKLT model has become a paradigm for the problem of fixing moduli in flux vacua. While 

one can raise issues about the analysis (i.e. in what sense the approximations arc truly system- 
atic), it provides a plausible scenario for fixing of moduli in supersymmetric and approximately 
super symmetric IIB compactifications. It also provides an illustration of the ways in which 
realistic landscape states are unlike the BP model and similar theories which have been used 
to consider eternal inflation and related issues. 

The basic ingredients in KKLT are a IIB theory compactified on a CY Manifold (Orien- 
tifold). Before turning on fluxes, the light fields of such a compactification consist of complex 
structure moduli and dilaton {zi^r), and Kahler moduli {pa)- Fluxes yield a superpotential 
W{z,t). The equations D^^W = 0,DtW = 0, in general, fix Zi,T. 

For large p, the low energy effective theory consists of the Kahler moduli. In the case of a 
single Kahler modulus, KKLT argue that the superpotential and the Kahler potential have the 
form: 

W = Wo + e-^P; K = -31n(p + pt) (21) 
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If there are many states, then, in some, KKLT argued that Wq will be small; in these, the 
effective action gives rise to an AdS, supersymmetric minimum with p ~ — ■^ln(Wo)- Super- 
symmetry, they argued, could be broken by antibranes; alternatively, dynamics of low energy 
gauge interactions might provide a breaking mechanism. 

The KKLT model makes especially clear the limitations of the continuous flux approx- 
imation. The features of KKLT depend crucially on the smallness of Wq, which relies on 
cancelations of fluxes. For nearby choices of flux, Wq is order one in fundamental units. So in 
this instance, we need to step away from the continuous flux approximation, and consider what 
the neighboring states look like. These states might then be expected to have broken super- 
symmetry, with p fixed, if at all, at some high energy scale, and cosmological constant (positive 
or negative) also very large. Moreover, in all of these states, the potential vanishes for very 
large p. So how transitions might proceed towards the KKLT vacua is not clear. Transitions 
from relatively nearby states, might be expected to end, typically, in big crunches or with the 
universe rolling - eternally - to large p. It could well be that anthropic considerations permit 
one to ignore these "bad ends", but without a larger theoretical framework, it seems like this 
will be a challenging story to sort through. 

For the /^-symmetric vacua that we consider shortly, the situation may often be different. 
Configurations starting out in a finite neighborhood may well end up in the preferred state a 
finite fraction of the time. 

4.4 The Polonyi Model 

Setting aside the problem of p, if we consider models which capture more of the features of the 
IIB theory, we encounter a richer and more subtle structure. Consider, first, a familiar example, 
the Polonyi model. We will think of the parameters of the superpotential of this model as 
surrogates for flux (we will shortly consider compactifications with a similar structure, where 
these parameters are determined by fluxes). The model consists of a singlet chiral field z, with 
superpotential: 

W = a + Pz; K = zz*. (22) 

As one varies the parameters, this has both supersymmetric and non-supersymmetric minima. 
In the supersymmetric minima, {W) is not a simple linear function of the fluxes. 

-a ± ^o? - 4/32 a ± ^o? - 4/32 
' = 2^ ' = 2 
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For large /3, there are no supersymmetric solutions. The non-analytie behavior as a function 
of the fluxes is not surprising. If a = f3 = 0, there is a moduli space of vacua. Turning on the 
fluxes lifts the degeneracy. This phenomenon is analogous to phenomena in ordinary degenerate 
perturbation theory, in which non-analyticity in an expansion parameter is typical. 

The non-super symmetric branch exhibits similar behaviors. We can rewrite the superpo- 
tential as 

W = ail + uz) (23) 

with u = l3/a; a then plays the role of n in our earlier discussions, and u of one of the "angles" 
introduced there. To simplify the algebra, we will assume that a and /? are real; we will also 
assume that the vev of z, {z) = x is real (this preserves a CP symmetry of the lagrangian). 
Then the potential is 

V{x) = a^e^' (u^{x^ - + 1) + 2u{x^ - 2x) + x^ - s) . (24) 

The equation, = then takes the form: 

u^ix^ + x^) + u(2x^ - x^ - 2) + x^ - 2x = 0. (25) 

Solving for x in terms of u is difficult, but solving for u in terms of x is simple. One root is the 
supersymmetric branch we discussed above: 

The second, non-supersymmetric branch has: 

- = —^- (27) 

We focus here on the non-supersymmetric branch. Substituting u into the potential, yields 
the value of the energy at the minimum of the potential for a given x: 

^2 



X' 



6 



x^ - 8x2 4 



(28) 



Note that the vacuum energy is positive for: 

x^ > 4 + 2\/3 x^ < 4 - 2\/3. (29) 

We need to check stability; a little bit of algebra yields:: 

(PV 1 



dx^ x^ 



(X [-x"" + 4x'' + 12x^] . (30) 
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Note that this is positive for sufficiently smah a;, and badly behaved for large x. So we 
have established the existence of a range of parameters for which there are stable stationary 
points where the vacuum energy is positive and supersymmetry is broken. 

For large and small u, we can solve the equations. In both limits, 

la 

u = -- x = -- (31) 
For small /5 (large x), the vacuum energy is approximately 

V = (5-'e^"l^\ (32) 

while for small a (small x) 

V = 4/3Va^ (33) 

We see that the neighborhood of the a = /? = point is surrounded by a sea of negative 
cosmological constant states with a and (3 comparable. 

Thinking of a and (3 as fluxes, we will take these parameters as large integer multiples of 
some basic unit. We will estimate tunneling rates, as in [15], by assuming that the tension is 
proportional to Ax in the transition, and the bounce action is proportional to 

Sb ~ T^/{/\Vf (34) 

Then an interesting question to ask is: starting in a configuration with a <^ (3, does the system 
flow towards the region of AdS states with a /?, or does it flow towards the zero energy states 
with a = 0. Estimating the bounce action as above, it is easy to see that 

>gb(Aa = 1) ^ ^ 

5f,(A/3 = 1) a ^ ' 

i.e. the transitions which decrease /3 are faster than those which decrease a. The couplings 
tend to equalize, and one is driven to the negative c.c. regime. 

While a simple, toy, model, the Polonyi model already captures certain features of the IIB 
flux vacua. In particular, the expectation values of moduli fields are complicated functions of 
the fluxes (both for super symmetric and non-super symmetric vacua). As a result, the tunneling 
behavior of the system is complex, with the system, for example, tending to avoid the zero c.c. 
region and heading rather to a region with negative cosmological constant. We will encounter 
behaviors of this type in a more intricate model in the next section. 
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4.5 Multiple-Field Polonyi Model 



In the landscape framework, we are interested in models with many fields, analogs of the field 
z. We can ask whether there are typically mainly dS or AdS vacua. A model which is simple 
to analyze has P fields and a Sp symmetry: 

W = aJ2zi + b. (36) 

i 

To make the analysis simple we will: 

1. Look for stationary points where the Sp symmetry, and CP, are preserved, i.e. Zi = x. 

2. Take the Kahler potential to be z*Zi, and ignore the terms in the potential. 

We will study the potential in the limit of large P. In this limit, even though the resulting 
equations axe cubic, we will be able to construct a solution in powers of 1/ Vp. The potential 
is: 

V = J2i(^ + ZjW)ia + ZjW*)-3W*W. (37) 
j 

so 

dV ^ dW* dW 

g^ = W{a + z^W*) + ^(a + z*W)z,-^ - 3W—. (38) 
At the would-be minimum, 

W = aPx + b. (39) 

So, substituting above, we obtain 

= P2(2a2x3) + P{3abx'^ - a^x) + {b'^x - 2ab) (40) 
For large P, in order to solve this equation, it is necessary that x take the form: 

+ ^ + . . . . (41) 



pl/2 ' p 

It is easy to see that " = As a result, at the minimum. 



W ^ ^PV2. (42) 
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On the other hand 



DW f^a + 



1 a 



pi/2 ^ 



(43) 



'2PV2 



a 



so the potential is: 




(44) 



So the solutions arc Dc Sitter, independent of a and b. Including e , one finds, for large P, 
that the minima lie near the origin, with cosmological constant 



Similar results hold for other choices of Kahler potential, such as = — ln(l + z*Zi). 

Clearly this is just a model, and one could well imagine that other parameter choices (and 
possibly other stationary points) would yield AdS minima, but this illustrates that there is no 
particular problem, with large numbers of fields, obtaining dS spaces. 

4.6 General Orientifolds of IIB on Calabi-Yau Spaces 

In order to relax some of the assumptions of the BP model, it is useful to consider flux vacua 
within the context of actual string compactifications in which the fixing of moduli is more or 
less understood. Perhaps the most widely studied flux vacua are provided by orientifolds of IIB 
theories on Calabi-Yau spaces. Here the candidate fluxes are RR and NS-NS three form fluxes. 
These are constrained by tadpole cancelation conditions. Most analyses have focussed on a 
low energy supersymmetric effective lagrangian, obtained by reduction of the ten dimensional 
supergravity theory. This lagrangian is described by a Kahler potential and superpotential for 
"light" fields. There are two classes of moduli which must be understood: complex structure 
moduli, Zi, and Kahler moduli, pa- Classically, the superpotential depends only on the Zi 
fields. Typically, one first tries to solve the cqiiations for a supersymmetric stationary point 
with respect to the complex structure moduli, Dz^W = 0. The Kahler moduli yield a no-scale 
structure, with vanishing potential. If the solutions of the equations Dz^W = yield W = 0, 
then the classical system is ftiUy supersymmetric. One might expect, for complicated models, 
that W = is the result of an (unbroken) discrete R symmetry. 

Our interest, for modeling transitions in the landscape, lies in non-supersymmetric solu- 
tions. There are two approaches one might adopt, closely related to analyses of Douglas and 
Denef [9] . In one, we will discard the Kahler moduli, keeping only the complex structure moduli. 



A Pa^. 



(45) 
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In this analysis, non-susy states arc stationary points of the potential, for which DiW 7^ 0. For 
most choices of flux, there arc, indeed, no supersymmetric solutions. These models resemble 
the Polonyi model, with many AdS minima. In the second, we keep the Kahlcr moduli, using 
the lowest order Kahlcr potential, and assuming that the Kahlcr moduli are somehow fixed at 
large values. SUSY is broken in all of these states, which are either dS or flat. 

4.7 Tq Inspired Model 

Perhaps the simplest example of this type is the oricntifold of the IIB theory on a lattice, 
studied by [22]. This theory is described by a set of complex structure moduli, Tij. Correspond- 
ingly, there are a set of NS-NS and R-R fluxes. In this section, we will consider a particular 
limit of this oricntifold, which will allow us to illustrate certain basic issues surrounding sym- 
metries. Our principle interest is in non-supersymmetric, dS vacua, which we will not be able 
to describe, in general, in any systematic way. So we will adopt a simple model. We will take 
the radius to be large and fixed, and examine the potential for the complex structure moduli, 
ignoring the dilaton. At the same time, we will include only the NS-NS fluxes. In this case the 
fields of the model are Tij, i,j = l,...3. They are described by a superpotential: 

W = ao det(T) -I- aij(cof (r))ij + bijTij + bo (46) 

and a Kahler potential 

K = ln{idet{T-T^)) (47) 

For special choices of the fiuxes, the resulting lagrangian exhibits global symmetries. For 
example, for 

ttij = ai5ij bij = biSij (48) 

the lagrangian possesses an SU{3) symmetry. This symmetry is broken by the underlying 
lattice to, at most, some discrete subgroup (we will discuss the symmetries of this model, in 
general, shortly), but this breaking is not visible in the classical theory. As we will see, such a 
choice of fiux is enough to fix all of the moduli. 

To summarize the model: 
1. We ignore the dilaton 
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2. We keep the Kahler moduli, pj, but assume that they are fixed at some large value, where 
they can be described by their classical Kahler potential, and W is independent of p. The 
effect of this is to cancel the 3|PFp in the lagrangian for the r^j's. 

3. We work in the continuous flux approximation, ignoring tadpole cancelation constraints. 

Even with these simplifications, the potential for the moduli is rather complicated to 
analyze. In certain limits, however, it is simple. An easy case to study is that with 60 = 0, and 
ai, bi proportional to the unit matrix. With r = ix, the potential is given by 

8x ^ ' 

which has a minimum for 



_ - 3aox^)(&i + aox^) 

At this minimum, one can check that there are no tachyons and the potential is positive provided 

-61 < aQX^ < ^^=. (51) 

17 + 015 

To understand the significance of this condition, it is useful to express in terms of ai, 61, ao- 

2 _ -{al + 2ao6i) ± ^(a? + 2ao6i)2 + I2albl 
- 6^2 • (52) 

4.8 Symmetries of the Model 

Before turning on fluxes, the Tg model has a large moduli space. On various subspaces, it 
exhibits symmetries. For example, for r proportional to the unit matrix, one has separate 
Z4 symmetries in each of three planes, as well as an S-^. In terms of the underlying (real) 
coordinates, x*, y*, i = 1, . . . , 3, the Z4's take 

^y' y'^ -x' (53) 

(i fixed). These are R symmetries, because they change the sign of the holomorphic three form, 

d^l = dz^ A dz^ A dz^ (54) 

with 

^ = 3^ + Tijy^ (55) 
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Odd permutations are also R symmetries. 

But we need to ask, as well, how the fluxes transform under the various symmetries. Ref. 
[22] provides a convenient catalog of three forms {H3{Te, Z)): 

ao = dx^ A dx'^ A dx^ Pq = dy^ A dy'^ A dy^ (56) 

aij = ^Eikidx^ A dx^ A dy^ ff^ = —^ejkedy'^ A dy^ A dx^ 
The corresponding fluxes are the quantities ao,bo,a^'' ,bij encountered above. 

It is easy to read off the transformation properties of the fluxes from these expressions under 
the various symmetries. Our particular choice, aij = aiSij, breaks all of the R symmetries. It 
leaves over cyclic transformations of the coordinates (i.e. even permutations). 

Now consider turning on small off-diagonal entries in oi, 61, breaking the microscopic cyclic 
symmetry (in the lagrangian itself, these break the SU{3) symmetry). For small ai, positive 61 
and negative ao, the condition, eqn. 51, for the absence of tachyons, reads, approximately: 

-h <-j< -bi/5 (57) 

which is always satisfied. Because there are no tachyons, this generates shifts in Tij of order 
Sai,Sbi, and changes in the energy of order 6a1,5b1. Transitions between states arc associated 
with five branes wrapping three cycles. If radii are of order one, then the tensions arc of order 
one (i.e. they are not enhanced by powers of flux). So, for small Aai, A61, transitions towards 
the symmetry point are slow; there is no particular tendency of the system to flow towards the 
symmetry point; only for very special initial conditions is symmetry achieved. 

4.9 R Symmetries within the Tq Model 

A simple example of an R symmetry in the Tg model is provided by taking a lattice which is a 
product of three square lattices in each of three planes. Consider the symmetry: 

xi -^yi yi^ -xi (58) 

This commutes with the orientifold symmetry. We can classify fluxes according to their behavior 
under this symmetry. The invariant three forms are: 

aia (59) 
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and 



aai (60) 

where a = 2,3. The fields ria, Tai transform hke the superpotential, i.e. by a phase e~ . With 
Til = i + 6t, St ^ —6t under the symmetry. 

To make our expressions simple, we will also impose a symmetry 2 3 on the fluxes. Then 
calling the fluxes (as in [22]) aia,bai, etc ,the superpotential is: 

W = ai2 [T21T33 — T23T31 + T31T22 — T32T2l\ + 021 [t12T33 — T13T32 + T13T22 " T23T12] + 

+ [ri2 + T13] + b^^ [t2i + T23] 

We are looking for supersymmetric stationary points with vanishing W, so we wish to solve 

dW 

(61) 



dTij 



with Til = Tai = Tia = 0. For general choices of the fluxes ai^, etc., the resulting equations, 
however, arc incompatible. In the language we used in section 3.3, there are more X type fields 
(ria,ra,i) than (p type fields (t22,T33). In more elaborate models (e.g. Calabi Yau's defined by 
intersections in weighted projective space) this is not the case. In the present model, solutions 
do exist for special choices of flux. In particular, if we take 

6^2 = 621 ai2 = a2i (62) 

then there is a moduli space with 

T22=T"33 T"23 = T"32 = T33 + — (63) 

ai2 

The next step is to turn on small symmetry violating fluxes, fixing these remaining moduli, 
and studying the behavior of the energy and fields as functions of the small fluxes. To first 
approximation, we simply need to write the effective action for the light moduli, and find its 
stationary points. There are many fluxes which can be turned on, and the resulting potentials, 
for which supergravity corrections are important, are complicated. To get some flavor for the 
problem, note that we need to fix rii,r22 and r33. So we can turn on small fluxes oq, and 
similarly 611,622,633 (again, to keep things simple, we will impose a 1 2 symmetry). Then 
the added terms in the superpotential are: 

Wi,e. = ao(riir22r33 - TUT23T32 + ...) + 6iirii + 622(r22 + T33). (64) 
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One can check that, with these choices, there are no tadpoles for ri^, a = 1,2. We can substitute 
for T23, T32 their forms on the moduh space, eqn. 63. The resulting expression is similar to that 
we studied in section 4.7, and we expect a similar structure for the solutions. Energies are of 
order (where the n's are the symmetry-breaking fluxes above). In particular, wc expect that 
the diagonal elements of r arc fixed, that for a range of fluxes there are no tachyons, and that 
transitions tend to keep the diagonal elements of a, b comparable. So shifts of the heavy fields 
will be of order v? jN'^ . 

Now we see that the tunneling rates are as we guessed for i?-symmetric stations in section 
3. Energies have the form: 

4 

£;o = ar? + 6-^ (65) 

so energy splittings between states of different TV behave as n^/N^. The tensions, again, are of 
order one, parametrically, the bounce actions for transitions which change N behave as 

S6(AAr = 1) ~ ^ (66) 
On the other hand, transitions which change n simply scale as 

Si,{5n = 1)7^\. (67) 
So symmetry-restoring transitions are highly favored. 



5 Conclusions 



Even if underlying fundamental physics is some type of landscape of metastable vacua, it is 
unlikely that, in the near future, we will have anything resembling a complete catalog of these 
states. Instead, any sort of understanding of the connection of microphysics and the phe- 
nomenon we observe will have to be based on robust statements about distinctive features of 
physics at much lower energy scales. In an earlier paper [15], we focussed on supersymmetry 
and the question of stability. In the present paper, we have mainly considered discrete symme- 
tries. We have remarked that states with discrete symmetries are rare[14], but suggested that 
they might be cosmological attr actors. This was motivated by the naive assumption that the 
potential might grow quadratically with flux around the symmetric point, so the system, start- 
ing in the neighborhood of the symmetric point, might tend to tunnel towards the symmetric 
point. To study this problem we have considered several model landscapes, which capture likely 
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features of a would-be underlying cosmic landscape. Even before approaching the problem of 
symmetries, we have argued that the Bousso-Polchinski model (the simplest model landscape 
we have considered) is likely too simple a picture, given its assumptions about positivity and 
constancy of flux coefficients, and given the existence of asymptotic pscudonioduli spaces. We 
have considered a number of toy models which capture more, but not all, of these expected 
features. 

To consider the question of whether symmetric points are attractors, we have worked in 
the continuous flux approximation. This is hardly likely to be realistic, but it should give some 
feeling of whether it is common for the system to reach symmetric points. Again, it provides 
a model, relatively simple to analyze, in which one can assess whether there is or is not a 
net drift towards the symmetric points. We have seen that plausibly an order one fraction of 
states lie within the neighborhood of symmetric points. However, we have also seen that for 
non-R symmetries, this fact is of no particular importance. Even if the system starts in the 
neighborhood, the system will drift towards non-symmetric, AdS states. 

For R symmetric states, the situation is potentially quite different. We have investigated, 
first, the question of whether there are likely to be predominantly dS or AdS states in the 
neighborhood. Within the models we have studied, both situations arise, so it seems plausible 
that, for a finite fraction of R symmetric points, the neighborhood is populated by metastable 
dS states. In this case, we have seen, that typical states in the neighborhood will drift towards 
a symmetric state. We take this to mean that an order one fraction of possible starting points 
leave the system in a symmetric state. Thus, even though such states are rare, cosmological 
considerations are likely to favor them. 

Based on this study, we offer the following tentative conclusions: 

1. An underlying landscape is likely to lead to low energy supcrsymmetry, simply because 
supersymmetry seems the most generic way to guarantee stability. 

2. Discrete R symmetries might plausibly emerge in a landscape. Even though these states 
are far less numerous than states without such symmetries, they are likely to be cosmo- 
logical attractors in a finite fraction of cases. As discussed in [19], they should tend to 
exhibit supersymmetry breaking at low energies, and there can easily be enough of them 
to resolve the cosmological constant problem. 

3. Discrete non-R symmetries (as one might want to understand, say, the structure of the 
NMSSM or some flavor structure) do not seem favored in any particular way. They are 
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rare, and don't seem to be singled out by (at least naive) cosmological considerations. 

4. Slightly increasing the size of the discrete symmetry does not seem to lead to a signifi- 
cant suppression, so intricate discrete symmetries (such as might suppress dimension five 
proton decay) might plausibly arise. 



As we have explained, establishing these statements as an ironclad "top down" prediction 
is difficult, but it seems likely that further landscape studies could well establish that they are 
robust. 

In this paper, wc have focussed principally on discrete symmetries. The problem of gauge 
symmetries is of obvious interest. In field theory, it is well-known that finite temperature tends 
to favor the maximal unbroken gauge symmetry. The issues, in the landscape context, are 
different and somewhat more complicated than those connected with discrete symmetry, and 
require a separate investigation, which will be reported elsewhere. 
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